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ABSTRACT 

We present results from 7 -ray observations of the Coma cluster incorporating 6 years of Fermi-LPfl data and 
the newly released “Pass 8 ” event-level analysis. Our analysis of the region reveals low-significance residual 
structures within the virial radius of the cluster that are too faint for a detailed investigation with the current 
data. Using a likelihood approach that is free of assumptions on the spectral shape we derive upper limits on 
the 7 -ray flux that is expected from energetic particle interactions in the cluster. We also consider a benchmark 
spatial and spectral template motivated by models in which the observed radio halo is mostly emission by sec¬ 
ondary electrons. In this case, the median expected and observed upper limits for the flux above 100 MeV are 
1.7 X 10“^ phcm“^ and 5.2 x 10“^ phcm“^ respectively (the latter corresponds to residual emission 
at the level of 1.8cr). These bounds are comparable to or higher than predicted levels of hadronic gamma-ray 
emission in cosmic-ray models with or without reacceleration of secondary electrons, although direct compar¬ 
isons are sensitive to assumptions regarding the origin and propagation mode of cosmic rays and magnetic field 
properties. The minimal expected 7 -ray flux from radio and star-forming galaxies within the Coma cluster is 
roughly an order of magnitude below the median sensitivity of our analysis. 

Subject headings: gamma rays: Coma cluster — gamma rays: cosmic ray interactions 
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1 . INTRODUCTION 

The radiative yields of energetic cosmic rays (CRs) travers¬ 
ing intracluster gas in galaxy clusters and interacting with 
background radiation fields span a wide spectral range from 
radio to high-energy 7 rays. Extended regions of radio emis¬ 
sion (referred to as halos and relics) have already been ob¬ 
served in many clusters (e.g., Ferrari et al. 2008). Comp¬ 
ton scattering of relativistic electrons by the Cosmic Mi¬ 
crowave Background (CMB) radiation is the dominant pro¬ 
cess for emission above 50 keV (where thermal emission from 
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hot intracluster gas is sufficiently weak) up to 0(100) MeV 
(e.g., Brunetti & Jones 2014). Searches for this non-thermal 
(NT) X-ray emission have not yielded conclusive results (e.g., 
Rephaeli et al. 2008; Ajello et al. 2009; Wik et al. 2009).“ 
At energies higher than 0(100) MeV the radiative decay of 
neutral pions (produced in energetic proton interactions with 
ambient protons, henceforth referred to as p-p) is expected to 
dominate, if CR interactions in the intracluster gas take place 
at non-negligible rates. 

While there is not yet observational evidence for energetic 
protons in clusters, the presence of relativistic electrons is 
well established from radio observations. Clearly, strong ra¬ 
dio galaxies in clusters are sources of relativistic electrons: 
e.g., M87 in Virgo (Bolton et al. 1949), NGC 4869 in Coma 
(Willson 1970), NGC 1275 in Perseus (van den Bergh 1961); 
see also Dutson et al. (2013) for a comprehensive search for 
the brightest cluster galaxies. It is also possible that ener¬ 
getic particles are (re)accelerated in the intracluster medium 
(ICM), by, e.g., merger and accretion shocks (for a recent re¬ 
view, see Brunetti & Jones 2014). From the large extent of 
cluster radio halos and typical synchrotron energy loss times, 
it is commonly thought that the emitting electrons are sec¬ 
ondary, namely decay products of charged pions (produced in 
p-p interactions; see, e.g., Brunetti & Jones 2014, for a recent 
review). Since contemporary viable models predict that the 
dominant 7 -ray emission process is 7 r° decay (e.g., Brunetti 
& Jones 2014) and since the electron lifetimes to radiative 
losses are generally significantly shorter than the source cross¬ 
ing time (Petrosian 2001), in our work here we do not consider 
purely leptonic models. 

An analysis of 50 clusters using four years of data from 
the Farge Area Telescope (FAT) onboard the Fermi satel¬ 
lite (Ackermann et al. 2014a) resulted in upper limits on the 
CR-induced 7 -ray emission (see also Ackermann et al. 2010; 
Huber et al. 2013; Griffin et al. 2014; Prokhorov & Chura- 
zov 2014). The 7 -ray upper limits from these analyses put 
stringent constraints on the energy density of energetic pro¬ 
tons in clusters, severely constraining models for acceleration 
by structure formation shocks (reviewed by Brunetti & Jones 
2014). 

Due to substantial interest in exploring cluster NT emission, 
and the availability of a longer FAT dataset with an improved 
event-level analysis (see below), a dedicated analysis of the 
most favorable cluster candidates is warranted. The Coma 
cluster is a natural choice for our search, being the nearest 
massive cluster with a bright radio halo. Coma is also located 
near the North Galactic pole where the diffuse gamma-ray in¬ 
tensity is at a minimum. The vast array of broadband obser¬ 
vations, including detailed measurements of the radio halo in 
Coma, provide a sound basis for testing theoretical models for 
the gamma-ray emission (e.g. Gabici & Blasi 2004; Reimer 
et al. 2004; Berrington & Dermer 2005; Pinzke &. Pfrom- 
mer 2010; Brunetti et al. 2012; Pinzke et al. 2015). Coma 
has been studied using the FAT (e.g., Ackermann et al. 2010; 
Arlen et al. 2012; Han et al. 2012; Prokhorov 2014; Zandanel 
& Ando 2014), its predecessor, EGRET (e.g., Reimer et al. 
2003), as well as Cherenkov telescopes (e.g., Aharonian et al. 
2009; Arlen et al. 2012), but has yet to be detected in MeV-to- 
TeV 7 rays. We report here on the deepest observation of the 
Coma cluster covering the MeV-GeV band to date, obtained 

Recent observations of the Coma and Bullet galaxy cluster with NuS- 
TAR only place upper limits on the intracluster component (Wik et al. 2014; 
Gastaldello et al. 2015). 
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using six years of LAT data analyzed with Pass 8 (Atwood 
et al. 2013). 

We describe the data analysis procedure in §2 and present 
our results in §3. A comparison of the results with predicted 
emission levels is presented in §4; we briefly summarize our 
conclusions in §5. In this paper we use the Planck mea¬ 
surement of the Hubble constant, Hq = 70kms“^ Mpc“^ 
(Planck Collaboration et al. 2014). At a distance of ~ 
100 Mpc (z = 0.023), the virial radius of the cluster of 2.0 
Mpc corresponds to a subtended angle on the sky, O 200 = 
1.23° (Ackermann et al. 2014a). 

2 . DATA SELECTION AND ANALYSIS 

The LAT is a pair conversion telescope with a large field 
of view of ^ 2.4 sr sensitive to 7 rays from ^ 20 MeV up to 
> 300 GeV (Atwood et al. 2009; Ackermann et al. 2012c).®^ 
This work uses 6 years (MET 239557414 - 428903014) of 
public Pass 8 LAT data. Pass 8 is an extensive rewrite of 
the core reconstruction algorithms (Atwood et al. 2013) ap¬ 
plied to the data taken by the LAT that is informed by in¬ 
flight performance and characteristics of the LAT, which re¬ 
sults in a increase in 7 -ray acceptance by 20 ^ 0 % with re¬ 
spect to source events reconstructed with Pass 7REP (Bregeon 
et al. 2013), depending on energy (see App. A for details). 
We select photons with energies from 100 MeV to 10 GeV 
within a region of interest (ROI) centered on the Coma clus¬ 
ter at 02000 = 194 . 95,(52000 = 27.98 and select a square of 
15° X 15° as our ROI.®^ 

Limiting the data selection to zenith angles less than 90° 
allows us to effectively remove photons originating from the 
Earth limb.®^ We use gtmktime to select time intervals dur¬ 
ing which the LAT was in nominal science operations mode 
and excluded intervals coincident with 7 -ray bursts and solar 
flares. We bin our data in 16 logarithmically spaced bins in 
energy and use a spatial binning of 0 ? 1 per pixel. 

In our analysis we include all sources listed in the lat¬ 
est Fermi catalog of point-like and extended sources (3LGL, 
Acero et al. 2015) along with the standard diffuse Galactic 
foreground emission model recommended for point source 
analysis. We include a spatially isotropic model accounting 
for the extragalactic diffuse 7 -ray background and misclassi- 
fied CRs. The normalizations of both Galactic and extragalac¬ 
tic diffuse emission models are left free in the likelihood fit. 
The model components are convolved with the parametrized 
detector response represented by the P8R2_SOURCE_V6 in¬ 
strument response functions (IRLs). 

We make use of recent developments of the Science Tools 
package that incorporate the finite energy dispersion of the 
LAT for the likelihood analysis by convolving the point 
sources in our ROI with the LAT energy dispersion.®® Note 
that the energy dispersion is already accounted for when cre¬ 
ating the above-mentioned diffuse models. During the like¬ 
lihood fit, we allow all sources that are separated from the 

Owing to the improvement with Pass 8 , the effective energy range has 
been extended both to lower and higher energies, allowing for a more efficient 
reconstruction of 7 rays up to a few TeV with respect to previous reconstruc¬ 
tion passes. At these energies however, the statistics remain small. 

^ The coordinates of the cluster center were taken from the NASA extra- 
galactic database (NED) https://ned.ipac.caltech.edu. 

Note that this choice is conservative to minimize the degeneracy of the 
extended emission from the cluster and the much brighter Limb emission. 

All resources discussed, including data, analysis software, source cata¬ 
log and diffuse models are available from the Fermi Science Support Center 
at http://fermi.gsfc.nasa.gov/ssc/ 


cluster center by less than 6 ? 5 to have a free normalization 
to allow for the broad PSL at the lowest energies (the 68 % 
containment radius of photons at normal incidence with an 
energy of 100 MeV is roughly 4°). This choice ensures that 
99.9% of the predicted 7 -ray counts (integrated over energy 
for an spectrum, which would roughly correspond to 

the predicted cluster emission) is contained within the chosen 
radius (that is substantially larger than (I 200 of the cluster, see 
§ 2.1 for details). 

In order to assess the completeness of our reference emis¬ 
sion model for the Coma region (without a cluster source), we 
use the gttsmap tool to search for any additional gamma- 
ray sources (Lig. 1). Lor each pixel in the map we evaluate 
the test statistic {TS), defined as twice the likelihood ratio 
between the best-fit value of the alternative hypothesis, in¬ 
cluding a source and the best-fit value of the null-hypothesis 
(background only). In the TS'-map a point source modeled as 
a power law with L = 2.3 is tested at each pixel of the map. 
We find four new point source candidates that coincide with 
regions with individual TS values > 25. The details of the 
improved background model along with the best-fit values of 
all free parameters are given in App. B. 

We find two residual structures whose peak TS values in 
individual pixels (with a pixel size of 0?05) are ~ 9 and ~ 13, 
respectively. These structures, located within the virial radius, 
are separated from each other by about 1 Mpc at the distance 
of the Coma cluster, if within the cluster. We further investi¬ 
gate this potential excess emission by repeating the TS-map 
calculation in a soft (E < 1 GeV) and a hard (E > 1 GeV) 
energy band. While the hard band reveals at least three dis¬ 
tinct areas of excess emission that spatially overlap with the 
coordinates of X Comae, and NGC 4839, respectively, the soft 
map does not permit an immediate association with discrete 
sources.®^ Moreover, when comparing our results with radio 
observations of the Coma cluster we find that the weak diffuse 
excess in the soft band roughly overlaps with that of the radio 
halo. In §4 we compare the level of this residual emission with 
predicted emission from likely sources in the cluster. We fur¬ 
ther investigate this apparent excess by creating spectral resid¬ 
uals in each energy bin z; for this we take the number of ob¬ 
served photons in each bin Ci and then evaluate the number of 
predicted photon counts given our background-only hypothe¬ 
sis irii and calculate the spectral residual = {ci — mf/rrii. 
We evaluate the residuals over the full ROI as well as within 
0.5 and 1.0 x 8200 show these results in Lig. 2. While the 
residuals evaluated over the entire ROI are practically negligi¬ 
ble, we note an indication of a slight excess within the cluster 
virial radius. 

2.1. Cluster Spatial Modelling 

In this section we motivate our benchmark spatial template, 
hereafter referred to as the cored profile, which is motivated 
by observations of the radio halo in Coma. 

Evidence for NT particle populations in clusters is cur¬ 
rently limited to measurements of extended regions of radio 
emission. The diffuse halo component in clusters (mostly 
intracluster emission) that remains after subtracting compact 
sources, has a power law spectrum with spectral flux f^, cx 
zz““, with a ~ 1.0 — 1.5 (Lerrari et al. 2008), considerably 
steeper than typical Galactic spectra (a ~ 0.60 — 0.75). The 
radio halo in Coma, which is centered on two dominant radio 

Note however that we conservatively do not attempt to model the emis¬ 
sion from these sources separately. 
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Figure 1. Significance map of the Coma region output from gttsmap {left: full energy range, middle: soft band, right: hard band). Each map has a dimension 
of X 4° and a resolution of 0?05 per pixel. The cyan-colored upright triangles denote positions of known blazars (Massaro et al. 2009) while the diamond 
and the square denote the position of X Comae and Coma A, respectively (with NED-based positions). The position of NGC 4839 is marked with an inverted 
triangle. The dashed circle corresponds to the angle subtended by the virial radius, 0200- The solid contours correspond to measurements of the Coma radio halo 
and relic using the Westerbork Synthesis Telescope (WSRT) at 352 MHz (Brown & Rudnick 2011). The WSRT observations were convolved with a 10' beam 
and have NVSS sources removed (although the residuals of the tailed radio source in the center still contribute significantly). The lowest contour corresponds to 
20 mjy beam” ^ (3cr) and increases in each contour in steps of 20 mJy beam” ^. 



10 ^ 10 ^ 10 "* 
Energy (MeV) 


Figure 2. Top\ Observed photons in 16 energy bins along with the total num¬ 
ber of predicted background model (comprising 3FGL sources together with 
extragalactic and Galactic diffuse emission templates) counts as a dashed 
(red) line. Bottom'. Spectral residuals determined from the full ROI, also eval¬ 
uated within the cluster virial radius (green and red markers, respectively). 
For visualization purposes, the curves corresponding to the residuals within 
the virial radius {R < 0.5 x r200 and R ^ ?’20o) are offset by 0.1 and 
0.2 dex in energy with respect to the black marker points. 


galaxies (NGC 4869 & NGC 4874), has a - 1.3 (Schlick- 
eiser et al. 1987; Thierbach et al. 2003). From the recent 352 
MHz WSRT measurements of Brown & Rudnick (2011) the 
halo is ^ 1 Mpc in size. The contour map in the latter pa¬ 
per, or, more conveniently, its conversion to an azimuthally 
averaged profile (as plotted in Fig. 10 of Planck Collaboration 
et al. 2013) can be parametrized by the following distribution 
in terms of the angular distance 9 from the cluster center 


fr[9) 




( 1 ) 


with 6 c — 15' and p ~ 1.4. 

The synchrotron spectral flux is proportional to the the pro¬ 


jected (along the line of sight) value of the product 
where Ne and B are the emitting electron (column) density 
and locally averaged value of the magnetic field, respectively; 
a is the spectral index of the radio-synchrotron emission. Due 
to heavy energy losses of primary electrons (with the required 
energy of more than a few GeV) propagating out of their 
source regions, it is usually assumed that the dominant con¬ 
tributions to the (observed) radio and (predicted) 7 -ray emis¬ 
sions are from charged (tt^) and neutral ( 7 r°) pion decays, 
respectively. The spectral intensity of the 7 -ray emission is 
predicted to have the characteristic 7 r° decay bump shape, un¬ 
less the emission is dominated by a Compton component due 
to scattering of the secondary electrons off the CMB, in which 
case a power law with an index nearly equal to a would be 
predicted. 

Whereas the spatial profile of Nc is essentially the same 
as that of both tt^ and 7 r°, respective sources of the radio 
and 7 -ray emission, the radio profile reflects also the spatial 
dependence of In the hot, fully-ionized intracluster 

gas, electrical conductivity is high, and magnetic fields are 
expected to be frozen into the plasma. If so, the field spatial 
dependence is expected to scale as with n being the intr¬ 
acluster gas density (Rephaeli 1988); the index in this scaling 
is only somewhat lower, 1 / 2 , if magnetic energy (rather than 
flux) is conserved. With a commonly used gas density profile 
of the form (1 + where Tc (oc 9c) is the gas core 

radius, and typically (3 = 2/3, it follows that the (projected) 
profile of the 7 -ray emission is expected to be considerably 
shallower than the measured radio profile.®^ 

Taking the radio profile from Eq. (1) and /3 = 2/3, the 
profile of the 7 -ray flux is predicted to be at most moderately 
steep, J^{9) oc (1 -I- 9"^ j9\)~^!. 

In addition to the cored profile, we consider models that 
bracket the two extremes of possible spatial models; one that 
is based on the cluster being modeled as a point source (point¬ 
like) and another in which we assume a uniform distribution 
of predicted photons out to the virial radius (disk-like). Note 
that the the tabulated data for the bin-by-bin likelihoods for 

For our purposes here, this profile is sufficiently close to the more realis- 
tic density distribution deduced in the analysis of Planck Sunyaev-Zel’dovich 
measurements of Coma (Planck Collaboration et al. 2013). 
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disk emission of varying size (0.1 x i ?200 — 1-0 x i? 20 o) is 
available in a supplementary tar package. 

2.2. Likelihood Analysis 

We use an extension to the standard LAT likelihood analy¬ 
sis similar to the calculation of the spectral energy distribution 
(SED) of a source, which we refer to as the bin-by-bin like¬ 
lihood. This approach allows to calculate flux upper limits 
in many narrow energy bins that can then be used to easily 
test various broadband theoretical models without the need 
for a dedicated likelihood analysis (Ackermann et al. 2014b). 
Because the sensitivity of the LAT to extended emission is 
expected to be substantially different than for point sources, 
we provide a set of three SEDs corresponding to different as¬ 
sumed spatial templates, following the considerations from 
§ 2 . 1 . 

Eor the spectral modeling, we assume the cluster emission 
in each energy bin to be characterized by a single power law 
with spectral index T = 2.3, and calculate the profile likeli¬ 
hood (Rolke et al. 2005) of the normalization parameter. Note 
that the values for the nuisance parameters are determined 
from a global fit over the entire energy range to avoid con¬ 
vergence issues. We make this fit prior to constructing the 
bin-wise likelihood. While the choice of the index in the bin- 
by-bin construction is somewhat arbitrary (even though moti¬ 
vated by that of the measured radio spectrum of the halo), it is 
found to have only a marginal (< 5%) effect on the resulting 
integral flux limits (as noted in the discussion of systematics 
in §3.1). 

3. RESULTS 

The resultant SED for the Coma cluster is shown in Eig. 3. 
The SED was deduced by first deriving the observed 95% con¬ 
fidence level one-sided upper limits, selected as the interval 
within which twice the difference in the log-likelihood with 
respect to the best-fit value of the alternative hypothesis (in¬ 
cluding the cluster) equals 2.71. These limits were then com¬ 
pared with those obtained by randomly selecting high Galac¬ 
tic latitude regions of the sky and repeating the analysis in 
these blank fields.®® The observed limits are typically within 
the 68% containment band, while a few bins fall in the 95% 
containment band. 

The information from each bin is combined into the global 
log-likelihood that contains the spectral information over the 
whole energy range of the analysis: 

= ( 2 ) 

i 

In this equation we have the binned Poisson likelihood term 
Ci, in each energy bin i, where Di corresponds to the ob¬ 
served counts in bin i and is the normalization (the only 
free parameter) of the power-law with T = 2.3 that we have 
assumed in each bin. 9 indicates that the nuisance parame¬ 
ters, 9, (i.e., normalizations of free point sources and diffuse 
background emission components) are taken from the global 
fit over the entire energy range. Irrespective of the origin of 

® We select 450 random positions in the sky with |fe| > 30° and exclude 
directions in which the center of the ROI coincides with either a 3FGL source 
or a cluster contained in the HIFLUCS catalog (Reiprich & Behringer 2002; 
Chen et al. 2007). For clusters as well as detected extended sources in 3FGL, 
we furthermore enforce that neither the center nor the innermost 4° are cho¬ 
sen when selecting a random sky position (see. Section VI in Ackermann 
et al. 2014b, for a detailed discussion on the subject of blank fields). 


the dominant emission process (be it hadronic or leptonic), 
in the relevant energy range explored here the spectrum can 
effectively be assumed to have a power-law form. Thus, we 
select a set of spectral indices for which we provide integral 
flux upper limits between 100 MeV and 10 GeV (bottom row 
of Eig. 3). The maximum TS value found in our analysis is 
~ 13 for the cored profile (T = 2.3). When using the blank 
fields to assess the null-hypothesis distribution, and consider¬ 
ing the trials factor associated with testing both for disk and 
point-like emission, this TS value corresponds to a global sig¬ 
nificance of ~ 1.8 ct. Given the low statistical significance of 
the residual emission, we determine integral flux upper limits. 

The deduced integral flux limits in the energy range from 
100 MeV to 100 GeV are 5.2 x 10“® phcm“^ for 
our benchmark model. Modeling the cluster as a point 
source yields a limit that is almost a factor 2 lower, 3.1 x 
10 “® phem”^ s“^, whereas for the uniform disk model we 
deduce a slightly higher limit 5.8 x 10“® phcm“^ s“^. Note, 
however, that the observed limit in each case is above the 
blank-field median expectation, falling instead into the 68% 
containment band for the disk and the 95% containment band 
for point-like and cored profiles, respectively. Our main re¬ 
sults are summarized in Table 1. 

The limits we present here appear to be weaker than previ¬ 
ously reported ones, most notably when comparing our results 
with those from the most recent analysis of the Coma cluster 
using Fermi-LPLA data (Zandanel & Ando 2014). However, 
the observed limits (solid black line in Eig. 3) cannot be com¬ 
pared naively this way because our analysis uses a larger (and 
different) dataset and an improved background model in terms 
of sources included in the ROI (see, App. A for details re¬ 
garding a comparison of five years of Pass 7 and Pass 8 data). 
Since Zandanel & Ando (2014) report limits with associated 
TS values of zero, these results should instead be compared to 
our expected sensitivity based on the blank field study. Doing 
so we find that the expected sensitivity for emission from a 
point-like source is similar to that reported in earlier works by 
Zandanel & Ando (2014), but the sensitivity for an extended 
source has improved significantly; the median expected limit 
for the disk-like emission for T = 2.0 is a factor two lower 
than what was reported in earlier works. 

3.1. Systematic Uncertainties 

To assess the robustness of our results we perform a num¬ 
ber of systematic checks and vary fiducial parameters such as 
the spectral or spatial binning as well as the number of free 
sources and repeat our analysis for these choices. In partic¬ 
ular we investigate how these changes affect the bin-by-bin 
likelihood evaluations and how these changes influence our 
likelihood analysis for soft and hard spectra. We tested our re¬ 
sults against a set of IREs that represent minimal or maximal 
boundaries in the computation of the effective area and PSP 
within the systematic uncertainties of Pass 8 that are chosen to 
maximize and minimize effective area and PSP, respectively. 
To explore the uncertainties associated with the diffuse fore¬ 
ground emission, we follow the method in Ackermann et al. 
(2014a) by employing a set of alternative models (Ackermann 
et al. 2012b; de Palma et al. 2013). The total systematic uncer¬ 
tainties of the analysis we present here typically do not exceed 
22% for E > 300 MeV. The fractional uncertainty increases 
up to 54% if the full energy range (100 MeV ^ ^ 10 GeV) 

is considered. Eor extended sources with fluxes near the cur¬ 
rent sensitivity limit of the LAT, the uncertainty from diffuse 
foreground modeling is at least as large as the uncertainties 
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Figure 3. Top row: Black arrows are the observed 95% C.L. one-sided upper limits on the integrated energy flux per energy bin when modeling the cluster as a 
point source left, and as an extended source by considering either the cored profile, with a core radius of 0.3° {middle), or a uniform (‘disk’) profile with radius 
0200 = 1-23° (right). In each of these plots, the green and yellow bands indicate the 68 % and 95% containment bands, respectively, obtained by repeating the 
analysis at 450 random high Galactic latitude locations in the sky. These bands provide an estimate of the overall analysis sensitivity given the assumed emission 
model. The dashed black curve in both panels indicates the median value. Bottom row: Observed integral flux limits over the entire energy range from 100 MeV 
to 10 GeV for various spectral indices assuming a single power law in 7 -ray flux for the three cluster models. The data for the full profile likelihood for each 
model and in each energy bin are available in the html journal. 


Table 1 

Summary of the Likelihood Analysis Results 


Spatial Model 

Power-law index 

> too MeV) 
(xlO-^phcm-^s-l) 

> too MeV) 
(xlO-Sphcm-^s-l) 

Cored Profile 

2.3 

5.2 

1.7 

Point Source 

2.3 

3.1 

1.1 

Disk 

2.3 

5.8 

2.1 

Cored Profile 

2.0 

3.2 

1.1 

Point Source 

2.0 

1.7 

0.5 

Disk 

2.0 

3.8 

1.4 


Note. — From left to right: spatial model assumed along with the adopted power-law index followed by the observed integral flux upper limit. The last column 
denotes the median expected integral flux upper limit that we determined from blank fields. All limits are calculated at 95% C.L. Note that while we refer to 
r = 2.3 as our benchmark model, we list also the values for F = 2.0 for easier compaiison with previous works (see text for details). 


from the other contributors, as shown in Tab. 1. 

4. DISCUSSION 

Our analysis of 6 years of LAT Pass 8 data above 100 MeV 
for the Coma cluster region does not reveal excess emission 
at a significance level that is large enough to claim detection 
of diffuse emission in the ICM. The TS map in Figure 1 is 
nonetheless of appreciable interest; there seems to be residual 
(background-subtracted) emission from an area that overlaps 
partly with the Coma virial radius. It is therefore reasonable to 
compare the sensitivity of our analysis (as gauged through the 
blank field analysis) with the emission predicted from likely 
sources in the cluster. 


The deduced median integral flux upper limits for our 
benchmark and disk models with T = 2.3, which corre¬ 
sponds to the measured index of the radio (spectral energy) 
flux specified above (a = 1.3), translate to the luminos¬ 
ity bounds of L{> 100 MeV) ~ 1.4 x 10^^ ergs“^ and 
L{> 100 MeV) 7 ^ 1.8 x lO'^^ ergs“^, respectively, on the 
combined emission from all compact and extended sources in 
the Coma region. 

The combined 408 MHz flux density from the two central 
radio galaxies NGC 4869 (5C 4.81) and NGC 4874 (5C 4.85) 
is 2.1 X 10“^^ ergcm“^ Hz“^ (Willson 1970; Kim et al. 
1994). An estimate of the combined 7 -ray luminosity of these 
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Table 2 

Budget of systematic uncertainties 


Category 

Variation 

SED impact 

Integral (hard) 

Integral (soft) 

IRFs 

bracketing MIN/MAX** 

(cTAefF ~ 3%, CTpSF ~ 5%) 

< 10 % {E < 300 MeV) 

< 1% (300 MeV <E< 2GeV) 

< 5% (E > 2GeV) 

< 3 % 

< 5 % 

Diffuse modeling 

alt. diffuse models*’ 

< 50% {E < 300 MeV) 

< 10%(t; > 300 MeV) 

< 12 % 

< 28% 

free radius‘s 

10 deg (nominal 6.5) 

< 10% 

< 4 % 

< 6% 

spectral binning 

±50% 


< 10 % 

< 10% 

spatial binning 

±50% 

< 10% 

< 2 % 

< 2% 

SED index 

±10% 

< 2% {E > 400MeV) 

< 10%(t; < 400MeV) 

< 5 % 

< 5% 

Total** 


< 54% {E < 300 MeV) 

< 22% (E > 300 MeV) 

< 21 % 

< 42% 


Note. — For the evaluation of systematic uncertainties we consider two cases: the effect on the bin-by-bin likelihood (here the uncertainties should be applied 
to the upper limit reported in each bin) and the effect on the broad-band spectra assuming a hard (F = 1.6) and soft (F = 2.6) power-law spectrum when all 
spectral bins are combined (here the uncertainty is given with respect to the integral flux upper limit). 

“ See § 5.7 and § 6.5 in Ackermann et al. (2012a) for a detailed discussion of the bracketing IRF approach 
See Ackermann et al. (2012b) for the detailed model description. 

Compared to the nominal choice of 6?5, there are 13 more sources that are left to vary in the fit. 

Values have been added in quadrature. 


two galaxies can be readily made if it is assumed that the ra¬ 
dio continuum emission is mostly due to Compton scattering 
of the radio-producing electrons by the CMB, a reasonable 
expectation in light of the conclusion reached by Abdo et al. 
(2010) in their analysis of Fermi-LAT measurements of Cen 
A, the nearby radio galaxy. If so, the total Compton lumi¬ 
nosity, Lc, can be estimated from the measured radio (bolo- 
metric synchrotron) luminosity, Ls, by the simple scaling 
Lc ^ LsPo/P b, where po and pB = are the (present) 

CMB and magnetic field energy densities, respectively. With 
a mean field strength of ~ 5 p,G across the the extended radio 
regions of NGC 4874 and NGC 4869 (Feretti & Giovannini 
1987; Feretti et al. 1990), we compute Lc ~ 6 x 10"^° erg s“^, 
and L(0.1—10 GeV) ~ 2x10'*° ergs“*. Additional emission 
from other radio galaxies in the cluster renders this estimate a 
lower limit on the total 7 -ray emission from all radio sources 
in Coma. 

The other discrete (‘compact’) sources are mostly from 
star-forming galaxies, whose contribution to the cluster 7 - 
ray emission was estimated by Storm et al. (2012) based on 
a scaling relation between galactic IR and 7 -ray emission. 
They estimated that the superposed 0.1 — 100 GeV luminos¬ 
ity of star-forming galaxies in Coma is in the wide range of 
~ 3 X 10^° — 3 X 10"*^ erg s“*. Together, the estimates of ra¬ 
dio and star-forming galaxies yield a lower limit on the com¬ 
bined emission from galaxies in Coma, Lgai(> 100 MeV) ~ 
5 X 10"*° erg s“*. (We note in passing that the predicted 7 -ray 
emission from X Comae, a background AGN at 2 ; = 0.091 lo¬ 
cated about a degree north to the Coma center, is too weak to 
be of relevance for our discussion here, since its radio flux is 
much weaker than the two radio galaxies within the cluster.) 

Quantitative estimates of the intracluster 7 -ray emission 
can only be made in the context of specific models for intr¬ 
acluster energetic protons and electrons producing this emis¬ 
sion in interactions with intracluster gas (via 7 r° decay) and 
Compton scattering off the CMB, respectively. A viable 
model must include the particle source distribution, the prop¬ 
agation mode, and the spatial distributions of intracluster gas 
and magnetic field. While models for cluster NT particles 


abound (reviewed recently by Brunetti & Jones 2014), and 
their predicted levels of 7 -ray emission span a wide range, 
we can roughly estimate a minimal level of 7 -ray emission 
from the measured radio halo luminosity, avoiding the need 
for a detailed model (and inherent untested assumptions). The 
halo flux is about twice (at 408 MHz) that of the two central 
radio galaxies, and the mean (volume-averaged) intracluster 
magnetic field is at least a factor five lower than the value 
adopted for the central radio galaxies (Bonafede et al. 2010). 
We would then expect that the total number of radio-emitting 
electrons in the halo to be considerably higher than that in the 
central radio galaxies. In the context of a secondary electron 
origin of the measured radio halo emission, it can be shown 
that the extended 7r°-decay 7 -ray emission is comparable to or 
higher than that from the central radio galaxies. Indeed, this 
has recently been quantified in a detailed study by Rephaeli & 
Sadeh (2016, in preparation). 

The above considerations lead to the lower limit L{> 
100 MeV) ~ 1 X 10** ergs * on the combined galactic and 
extended 7 -ray emission from Coma. Based on this estimate 
the predicted minimal 7 -ray flux is at a level that is an order 
of magnitude lower than our deduced upper limit. Aside from 
the instrument response functions, the measurement sensitiv¬ 
ity is mainly governed by two factors: the detected (and thus 
modeled) discrete sources observed by Fermi, i.e., the 3FGL, 
and the model for the Galactic and isotropic 7 -ray foreground 
emission. Moreover, the residuals found in this work may be 
the first indication of cluster emission which could become 
more significant with further observations. 

The extrapolated sensitivity (given by the median expecta¬ 
tion from blank fields, scaled to a longer exposure time) of our 
analysis to 10 years of LAT exposure is comparable to recent 
predictions of diffusion models of CR protons with turbulent 
reacceleration of leptonic secondaries (Brunetti & Lazarian 
2011; Brunetti et al. 2012). However, even if the exposure of 
the object could be doubled by the end of the Fermi mission, at 
least in the context of the analysis techniques presented here, 
it appears unlikely that Coma will be significantly detected in 
the LAT data given the current results. This forecast, based 
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upon the enhanced sensitivity of Pass 8 and a comparison to 
blank fields, is in agreement with the conclusions of Zandanel 
& Ando (2014), who found no indications of residual emis¬ 
sion (TS' ~ 0) in an analysis of 5 years of Pass 7 reprocessed 
LAT data. 

5. SUMMARY AND CONCLUSIONS 

The main conclusions from this paper can be summarized 
as follows: 

1. In an analysis of 7 rays between 100 MeV-10 GeV col¬ 
lected over a period of six years with the LAT and re¬ 
processed with Pass 8 , we find excess emission within 
the cluster virial radius; however, the statistical signif¬ 
icance of this emission is well below the threshold to 
claim detection of 7 -ray emission from the cluster. 

When spectrally analyzed, this excess emission above 
the background expectation can be separated into a soft 
and a hard component; the soft {E < \ GeV) ap¬ 
pears to roughly spatially overlap with parts of the well 
measured giant radio halo, while the hard component 
(E > \ GeV) can be associated with parts of the halo. 

2. Using a 7 -ray template derived from the combination 
of Planck and WSRT observations of the radio halo, we 
find a maximum TS value ~ 13 for power-law emis¬ 
sion with a spectral index ~ 2.3 (after correcting for 
trial factors, this corresponds to a global significance of 
~ 1.8 cr). We derive limits on the integral CR-induced 
7 -ray flux, and our observed limit excludes fluxes above 
5.2 X 10“® phcm“^ s“^. While we focus in this paper 
on hadronic models where the radio emission is of sec¬ 
ondary origin, we emphasize that the results presented 
here are of a more universal nature; the relatively small 
variation of ~ 60% between our benchmark model and 
the most extreme disk profile indicates that our results 
are robust with respect to the assumed spatial distribu¬ 
tion of 7 rays. Our results and the bin-by-bin likelihood 
profiles provide the basis for comparisons with a vari¬ 
ety of specific models that will help improve our under¬ 
standing of CR physics in the Coma cluster.^*^ 

3. Based on scaling considerations and radio measure¬ 
ments, we derive a robust lower limit on the 7 -ray lumi¬ 
nosity of Coma that is a factor ^ few below the median 
sensitivity given our analysis of blank fields. 
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within the Coma ROI of 15° X 15°. Note that for Pass 7 we apply a different zenith angle cut (100°) but otherwise the basic ROI selection cuts are the same. 
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APPENDIX 

A. PASS 8 IMPROVEMENTS 

Pass 8 is the successor of the previous event level analysis (Pass 7REP) (Ackermann et al. 2012a; Bregeon et al. 2013). The 
improvements include a tree-based pattern recognition algorithm to identify and reconstruct tracks, a new energy reconstruction 
that uses minimum spanning trees and calorimeter clustering, which also extend the energy range at both high and low energies 
and an improved background rejection algorithm from the anti-coincidence shield. Pass 8 also improved on the implementation 
and training of the classification trees used in the gamma-ray selection algorithms, resulting in increased gamma-ray selection 
efficiency with respect to Pass 7, while keeping the same background rejection power (Atwood et al. 2013).^^ Eor high-level 
science analysis this presents a significant improvement in all metrics. Among the key improvements are a ^ 25% increase in 
the acceptance (and more than 50% below 100 MeV and above 300 GeV) along with an improved angular resolution. As a result 
the point source sensitivity in the energy range between 1 GeV and 10 GeV is enhanced by 30-40% with respect to Pass 7REP. 
Similarly, the improved PSP provides increased sensitivity towards spatially extended sources. 

Due to the nature of the changes introduced with Pass 8 , there is only partial overlap between the Pass 7REP and Pass 8 event 
samples, with the Pass 8 sample generally being the larger of the two. The reason for this is two-fold; the changes in the event 
selection cause different events to pass the source selection criteria in Pass 8 than in Pass 7REP as well as selecting different 
residual non-photon events. The increase in 7 -ray acceptance on the other hand increases the number of reconstructed events. 
Pig. 4 shows the fraction of shared events between our here reported Pass 8 dataset and the older Pass 7REP event sample after 
subjecting both datasets to the same ROTbased selection criteria. The total number of photons in our dataset is about 612k. 
Below ~ 3 GeV, the fraction of shared events drops considerably, and towards the lowest energies there are about four times 
as many events in the Pass 8 dataset that can partially account for the differences between previously published results based on 
Pass 7 and the analysis we present here. 


B. DETAILS OP BACKGROUND MODEL 

The basis for our background model is the 3PGL catalog of sources along with a set of diffuse templates rescaled to be used with 
Pass 8 . We find several grid positions in the significance map (described in Section 2) that correspond to statistically significant 
(TS > 25) 7 -ray excesses over the initial background model. The coordinates of four new point sources added to the background 


Performance plots for Pass 8 and its 
comparison with Pass 7REP are provided at 


http://www.slac.stanford.edu/exp/glast/groups/canda/lat_Performance.htm. 
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Table 3 

Point Source Candidates in ROI 


Name 

R.A. 

Dec. 

Model 

Parameters 


TS 


(°) 

(°) 





xFGL J1904.00-H3465 

190.40 

34.65 

PL 

r = 1.88, Eo = 3496 MeV 

112 

xFGL J1927.34-H3129 

192.73 

31.29 

LP 

o = 1.83, /3 = 0.39, Eh = 

2407 MeV 

86 

xFGL J2027.04-H2955 

202.70 

29.55 

LP 

a = 2.57, B = 0.41, Et = 

1237 MeV 

280 

xFGL J1914.49-H2080 

191.45 

20.80 

PL 

r = 3.45, Eo = 1000 MeV 

1220 


Note. — From left to right: name, right ascension, declination and assumed model along with fixed 
spectral parameters and TS values obtained in likelihood fit to the Coma ROI. All coordinates are given 
in J2000 epoch. 


Table 4 

Best-Fit Parameters of Background Model 


Source Name Normalization, No 

(cm~2 MeV“^) 


3FGLJ1230.3-H2519 

3FGLJ1231.7-H2847 

3FGLJ1254.5-H2210 

3FGLJ1258.1-H3233 

3FGLJ1258.4-H2123 

3FGLJ1301.5-H3333 

3FGLJ1303.0-H2435 

3FGLJ1310.6-H2446 

3FGLJ1310.6-H3222 

3FGLJ1314.8-H2349 

3FGLJ1321.0-H2215 

3FGLJ1323.0-H2942 

3FGLJ1326.1-H2931 

3FGLJ1332.8-H2723 

XFGLJ1914.49-I-2080 

xFGL J1904.00-I-3465 

XFGLJ1927.34-I-3129 

xFGL J2027.04-I-2955 

Extragalactic Diffuse” 

Galactic Diffuse” 


(2.11±0.09) X10-12 
(0.93±0.04) X10-12 
(l.lOitO.lO) X10-13 
(1.40±0.10) X10-12 
(0.30±0.10) X10-12 
(0.29±0.05) X10-12 

(2.80±0.10) X10-12 

(0.35±0.09) X10-13 
(2.60±0.05) X10-11 
(0.53±0.03) X10-12 
(4.20±0.10) X10-12 
(0.60±0.03) X10-12 
(0.40±0.30) X10-11 
(0.66±0.09) X10-12 
(1.50±0.20) X10-13 
(1.40±0.40) X10-11 
(O.eOitO.lO) X10-13 
(2.10±0.40) X10-13 
(1.04±0.01) 
(1.04±0.02) 


” The fitted value corresponds to the overall (unit¬ 
less) normalization of an all-sky template. The 
nominal value is 1.0. 


model are given in Table 3, along with their parameters as derived from a broad-band spectral fit. For two sources we find that a 

LogParabola (LP) of the form ^ ~ ) gives a better fit to the data than a power law (PL); the remainder 

is modeled using a simple PL instead. The normalizations of each source is left free to vary in the likelihood fit. 

Table 4 shows the best-fit values for the background-only model, including our new candidate sources. The associated uncer¬ 
tainty refers to the 68% parabolic error reported by MIGRAD. Note that the background model also comprises sources that are 
outside the selected 15° x 15° ROI, such as the bright source 3FGL J1224.9-I-2122 (4C -1-21.35), whose parameters are fixed to 
the reported values in 3FGL. 



